Quantum gravity corrections to neutrino propagation 
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Massive spin- 1/2 fields are studied in the framework of loop quantum gravity by considering a 
state approximating, at a length scale C much greater than Planck length £p = 1.2 x 10~ 33 cm, a 
spin-1/2 field in flat spacetime. The discrete structure of spacetime at £p yields corrections to the 
field propagation at scale C. Next, Neutrino Bursts (p « 10 5 GeV) accompaning Gamma Ray Bursts 
that have travelled cosmological distances, L « 10 10 l.y., are considered. The dominant correction is 
helicity independent and leads to a time delay w.r.t. the speed of light, c, of order (p £p)L/c « 10 4 s. 
To next order in p £p the correction has the form of the Gambini and Pullin effect for photons. Its 
contribution to time delay is comparable to that caused by the mass term. Finally, a dependence 
Log 1 oc p 2 £p is found for a two-flavour neutrino oscillation length. 
£ ■ PACS: 04.60.Ds 14.60.Pq 96.40.Tv 98.70.Rz 

The fact that some Gamma Ray Bursts (GRB) originate at cosmological distances, (~ 10 10 light years) 0, together 
with time resolutions down to submillisecond scale achieved in recent GRB data , suggest that it is possible to probe 
fundamental laws of physics at energy scales near to Planck energy Ep = 1.2 x 10 GeV ||f|- Furthermore, sensitivity 
will be improved with HEGRA and Whipple air Cerenkov telescopes and by AMS and GLAST spatial experiments. 
Thus, quantum gravity effects could be at the edge of observability |3|,[|]. Now, quantum gravity theories imply 
different spacetime structures and it can be expected that what we consider flat spacetime, can actually involve 
■ dispersive effects arising from Planck scale lengths. Such tiny effects might become observable upon accumulation 
| over travels through cosmological distances by energetically enough particles like cosmological GRB photons. 

Now, the most widely accepted model of GRB, so called fireball model, predicts also the generation of 10 14 — 10 19 
eV Neutrino Bursts (NB) [pra. Yet, another GRB model based on cosmic strings requires neutrino production M. 
Present experiments to observe high energy astrophysical neutrinos like AMANDA, NESTOR, Baikal, ANTARES and 
Superkamiokande, for example, will detect at best only one or two neutrinos in coincidence with GRB's per year. The 
planned Neutrino Burster Experiment (NuBE) will measure the flux of ultra high energy neutrinos (> 10 TeV) over 
^ | a ~ 1km 2 effective area, in coincidence with satellite measured GRB's |Tc| ]. It is expected to detect s=s 20 events per 
I | year, according to the fireball model. Hence, one can study quantum gravity effects on astrophysical neutrinos that 
1 ■ might be observed or, the other way around, such observations could be used to restrict quantum gravity theories. 

In this letter, the loop quantum gravity framework is adopted. In this context, Gambini an Pullin studied light 
propagation semiclassically |^|. They found, besides departures from perfect non-dispersiveness of ordinary vacuum, 
helicity depending corrections for propagating waves. In the present work, the case of massive spin-1/2 particles 
in loop quantum gravity is studied also semiclassically. They could be identified with the neutrinos that would be 
produced in GRB. Central ideas and results are presented whereas details will appear elsewhere [ pr[ . 

Loop quantum gravity H uses a spin networks basis, labelled by graphs embedded in a three dimensional space 
E. Physical predictions hereby obtained are a "polymer- like" structure of space jl2| and a possible explanation of 
black hole entropy A first attempt to couple spin--| fields to gravity, along these lines, was made in |Q and a 
generalization to the spin networks basis has been developed in |l5| . A significant progress in the loop approach to 
quantum gravity was made by Thiemann, who put forward a consistent regularization procedure to properly define 
the quantum Hamiltonian constraint of the full theory, which includes the Einstein plus matter (leptons, quarks, Higgs 
particles) contributions [jl6). It is based on a triangulation of space with tetrahedra whose sides are of the order ip. 
The cornerstone of Thiemann's proposal is the incorporation of the volume operator as a convenient regulator, since 
its action upon states is finite. Having at our disposal a regularized version of the quantum Hamiltonian describing 
fermions coupled to Einstein gravity, we will further need a loop state which approximates a flat 3-mctric on E, at 
scales £ much larger than the Planck length. For pure gravity this state is called weave p7]]. A flat weave \W) 



*jalfaro@chopin.fis.puc.cl 

t Associate member of Abdus-Salam ICTP, Trieste, Italy, hugo@xanum.uam.mx 
*me@aurora.nuclecu.unam.mx 



1 



is characterized by a length scale C >£> lp, such that for distances d > C the continuous flat classical geometry is 
regained, while for distances d <K C the quantum loop structure of space is manifest. The stronger the unequality 
d » lp holds, the more isotropic and homogeneous the weave looks like. For example, the metric operator g a b 
satisfies (W^gablVF) = 5 a b + O(^f-). Now, a generalization of such an idea to include matter fields is required. For 
our analysis it suffices to exploit the main features that a flat weave with fermions must have: in particular, it must 
reproduce Dirac equation in flat spacetime and this is just the basis of our approximation scheme, ft is denoted by 
|W, £), has a characteristic length L and it is referred to simply as a weave. 

The use of Thiemann's regularization for Einstein Dirac theory naturally allows the semiclassical treatment here 
pursued; expectation values w.r.t. \W,£) are considered thereby. They are expanded around relevant vertices of the 
triangulation and a systematic approximation is given involving the scales lp <5C Ad <SC Ac, the last two corresponding 
to, respectively, De Broglie and Compton wavelengths of a light fermion. Corrections come out at this level. 

The Hamiltonian constraint for a spin-^ field coupled to gravity consists of a pure gravity contribution, a kinetic 
fermion term, namely, 

ff« := f d s x ff (m T nV a i + c.c.) , (f ) 

and other terms Jl6[ whose contribution is summarized in (Q) below. Their analysis is an extension of the one for 
given here and will be spelled out in We use f = — ^a, the latter being the standard Pauli matrices. 

The fermion field is a Grassmann valued Majorana spinor ^ T — (ip T , (— ia 2 ip*) T ). The two component spinor tp has 
definite chirality and it is a scalar under general coordinate transformations. Hence, (0) is not parity invariant. The 
configuration variable is £ = {det{q))^ 1 which is a half density. The corresponding momentum is, with this choice, 
7r = similarly as in flat space. The gravitational canonical pair consists of Ei a and the su(2)-connection A J b (T>), 
where (detg)g ab = Ei a E ib . 

Upon regularization |T^| , the expectation value of (0) with respect to the weave becomes 



Wf|A<"M = -4 E ^* E ■ 
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(W, ClCs (v + s K (A)) ^A-y (r k h SK (A) ) ab wua (v)w jJA (v)\W, 



d_ 

■dj(v)> 



-(W,£\T k av) 1 ^\w iIA (v)w J j A (v)\W,0 -c.c. \. (2) 



Here an adapted triangulation of E to the graph 7 of the weave state |W, £) is adopted. Auxiliary quantities used 

are WkiA = Tr (rkh Sl ( A ) ^/(A)' V^Ki" ^ , where V v is the volume operator restricted to act upon vertex v. /i s (A) are 

holonomies along segments, s, of edges forming tetrahedra in the triangulation A Jig . ^(7) stands for the set of 
vertices of 7. The second sum, ^2 Sins jnsK=v , involves triples of segments s/, sj, sk intersecting at v. Notice that one 
actually averages over E(v) = n v (n v — l)(n v — 2)/6 possible triangulations (one for each triple of edges) when the 
vertex v is reached by n v edges (the valence) of the graph. 

To estimate (||) we associate to it c-number quantities respecting the index structure, together with appropriate 
scale factors arising from dimensional reasons and, most important, in line with the weave state approximating flat 
space with fermions. This amounts to an expansion of expectation values around vertices of the weave. The explicit 
form is taken from the expansion the involved operators would have in powers of the segments s a , \s a \ ~ £p -a 
procedure justified for weave states. Useful quantities coming in by expanding wua(v) = sj a Wi a + S] a sj b Wi a b + ■ ■ ., 
for instance, are 

W ia = \[Aia, a/KL ^iab = \tikl [Aka, [Ai b , \fV~v]] , (3) 

whose contribution to the average in the weave is estimated by considering that of Ai a and \/Vv to be of order ~ 1/C 

3/2 

and ~ lp , respectively. To proceed with the approximation we think of space as made up of boxes of volume 
£ 3 , whose center is denoted by x. Each box contains a large number of vertices of the weave, but is considered 
infinitesimal in the scale where the space can be regarded as continuous, so that we take C? « d 3 x. Let F(v) be a 

fermionic operator which produces the slowly varying (inside the box) funtion F(x); i.e. C « \d- Also let -A-G(v) be a 

'-p 

gravitational operator with average within the box G(x). The weave is such that YlveVij) E^y(^' £\F(. V )G( V )\W : — 
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Ebox(x) E«eB OX (S) 4 E^y(W.Cl4^ w )l W '^ = kd 3 xF(x)G(x). Notice that the tensorial and Lie-algebra 

structure should come out from flat spacetime quantities exclusively, i.e. °i? ia , r k , d b , e cde , g feim ; where 5 ah = °E ta °Ei b . 
In order to regain the flat spacetime kinetic term of the fermion Hamiltonian, we demand |W,£) to fulfill 



(W,Z\t B {v) 



d 



d£ A (v) 



(T k VP) AB w ta (v)w jb (v)\W,0 « (U B (v)Tr A (v)£p£ 



l f 2 {r k d^) AB °E ia {v) % b (v) 



(4) 



The second parenthesis here dictates the overall structure: indicates that each w ia (v) contributes a factor tp ^ , 

since the connection scales with 1/C (large length limit => flat spacetime), and W contributes a factor £p 3 ^ 2 . 
Independence on C of the final form of (0) gives the structure of the first parenthesis. The notation ^ stands for 
acting only upon £. By expanding (|^) at different orders in powers of s and using (|4|), one can systematically determine 
all possible contributions. Some examples of correction terms are 
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Sl a Sl b S j c SK d ir A (v)d d £b (v) (W, H(T k ) AB {w lab , w jc } \W, 



K5 £ 
i_C?_ 

A I 3 



d 3 x^{x)r k E™d d ax) 



E 



f 



e ijk e IJK a 8K c 8 d q ^q ^ { W ^\w id {v)w je {v)\W^) 

p — > b{v) 3! 



-ins Ip J d 3 x n(x)T k °E kc d c V 2 f (f) 
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A similar treatment can be performed for every contribution to the Einstein-Dirac Hamiltonian constraint 1 1 1 . It 
is important to stress that the prediction of the values of the corresponding coefficients Ki would require a precise 
definition of the flat weave (or even better a Friedman-Lemaitre- Robertson- Walker weave), together with a detailed 
calculation of the matrix elements. Instead, within the present approach, the neutrino equation, up to order £p, 
becomes 



d C 
ih—-ihAa-V + — 
at 2C 



A 



I + k v 



C 
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£(*, x)+m(a- (3 ih a • V) i a 2 C (*, x) = 0, 

£p 



f 4 v 2 



1 + K 8 - 
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Notice that K4 would produce an additional Dirac mass for the neutrino. Since we are considering particles with a 
Majorana mass m, we take K4 = 0. In contrast to |g|, we have found no additional parity violation arising from the 
structure of the weave. The dispersion relation corresponding to (m) is 



2a/3m :; 



(8) 



where A, B 7 C have been expressed in momentum space and depend on C. The ± in Eq. (|^) stand for the two neutrino 
helicities. Let us emphasize that the solution £(t, x) to Eq.(Q) is given by an appropriate linear combination of plane 
waves and helicity eigenstates, given that the neutrinos considered are massive. 

Typically, for neutrinos, \p> « Ac and our approximation is meaningful only if C < \p>- In this way we make sure 
that Eq.(^) is defined in a continous flat spacetime. From here on: K = c = 1. To estimate the corrections let us 
consider a massive neutrino with momentum p = p. A lower bound for them is obtained by taking 1/C ~ 1/ A^ = 
\p\ = p. Up to leading order in £p, we get 



E± (p) :=E±(p,C) \ p=Pi c=i/p~P- 



2p 



l P ((9 2 ± B A )f + (6 1 ± 3 )m 2 ) + (6 5 ± 9 e )£ 2 P p 3 , 



(9) 



where we are assuming that all 0i are numerical quantities of order one. Besides, these are known functions of the 
original parameters Ki plj . To leading order in p, the velocities are 
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v±(p) = J p L \ P = P< c=i/p = 1 - + Ki(Ip p) T ^-^{^- (10) 

The order of magnitude of the corrections arising from the present analysis is calculated using the following values: 
m = 1CP 9 GeV, p ~ 10 5 GeV, L = 10 10 l.y. = 0.5 x 10 42 (distance travelled by the neutrino from emission 
to detection on Earth). Let us observe that in Eq.(^), the ratio between second and first order contributions in 
lp behaves like (pip) ~ 10~ 14 . Now consider the gravitationally induced time delay of neutrinos travelling at 
velocities v± with respect to those traveling at speed of light: At v = \L(1 — v±)\ = \n%\L (p lp). Notice that this 
expression, though helicity independent, is of the same form as the one in Ref. for photons. In our case we 
obtain At v = 0.3 \k\\ x 10 4 s. Besides, this correction dominates over the delay due to the mass term ^ which is 
ps 10 _10 s. The second interesting parameter is the time delay of arrival for two neutrinos having different hclicitics: 
At± = L\(v+ —v~)\ = \K7\L (p lp) 2 pa 1.5|/«r| X 10~ n s. This correction is supressed by a factor of (p lp) with respect 
to the former and it is comparable to the time delay caused by the mass term. Finally, consider the characteristic 
length L os corresponding to two-flavour neutrino oscillations, given by L os — ^ E 2 l Eb ^ = where E a .b denotes de 
energy corresponding to the mass eigenstates of the neutrinos with masses m a ,b respectively. As usual, we assume 
that neutrinos are highly relativistic (p a » m a ) and also that p a ~ pf, = p ~ E. The phase $ s describing the 
oscillation is <& os = j^- , where L is the distance travelled by the neutrino between emission and detection. The energy 
difference for the corresponding two flavours is 

AE = + A Pl p 2 l P + A( P2 m 2 )l P ps (l0^ 26 + A Pl x 10~ 9 + 10" 40 ) GeV. (11) 
2p 

This result could yield bounds upon Ap x , which measures a violation of universality in the gravitational coupling for 
different neutrino flavours. For the above estimation, Am 2 ps 10 -21 GeV 2 , and A(p2in 2 ) ps Am 2 were used. Here pi 
are flavour dependent quantities of order one. To conclude, we notice that (|ll|) implies L^ 1 oc p 2 lp, seemingly an 
effect not considered previously 
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